arXiv:1503.00070v2 [astro-ph.GA] 1 Apr 2015 


Cloud-cloud collision as a trigger of the high-mass star formation; 

a molecular line study in RCW 120 

K. Torii 1 ’ 2 , K. Hasegawa 1 , Y. Hattori 1 , H. Sano 1 , A. Ohama 1 , H. Yamamoto 1 , K. 
Tachihara 1 , S. Soga 1 , S. Shimizu 1 , T. Okuda 3 , N. Mizuno 3 , T. Onishi 4 , A. Mizuno 5 , and Y. 

Fukui 1 

1 Department of Physics, Nagoya University, Chikusa-ku, Nagoya, Aichi 464-8601, Japan 
2 Sub-Department of Astrophysics, University of Oxford, Denys Wilkinson Building, Keble 

Road, Oxford 0X1 3RH, UK 

3 National Astronomical Observatory of Japan, Mitaka, Tokyo 181-8588, Japan 
department of Astrophysics, Graduate School of Science, Osaka Prefecture University, 1-1 
Gakuen-cho, Nakaku, Sakai, Osaka 599-8531, Japan 
5 Solar-Terrestrial Environment Laboratory, Nagoya University, Chikusa-ku, Nagoya 

464-8601, Japan 

torii@a.phys.nagoya-u.ac.jp 

Received_; accepted_ 




2 


ABSTRACT 


RCW120 is a Galactic Hu region having a beautiful ring shape bright in 
infrared. Our new CO J=l-0 and J= 3-2 observations performed with the 
NANTEN2, Mopra, and ASTE telescopes have revealed that two molecular 
clouds with a velocity separation of 20kins -1 are both physically associated with 
RCW 120. The cloud at —8kms _1 apparently traces the infrared ring, while the 
other cloud at —28kms _1 is distributed just outside the opening of the infrared 
ring, interacting with the Hll region as supported by high kinetic temperature of 
the molecular gas and by the complementary distribution with the ionized gas. 
A spherically expanding shell driven by the Hll region is usually discussed as the 
origin of the observed ring structure in RCW 120. Our observations, however, 
indicate no evidence of the expanding motion in the velocity space, being incon¬ 


sistent with the expanding shell moc 
by applying the model introduced by 


eh We here postula te an alternative that, 


Habe fe Ohtal (119921) . the exciting O star in 


RCW 120 was formed by a collision between the present two clouds at a colliding 
velocity ~30kms _1 . In the model, the observed infrared ring can be interpreted 
as the cavity created in the larger cloud by the collision, whose inner surface is 
illuminated by the strong UV radiation after the birth of the O star. We discuss 
that the present cloud-cloud collision scenario explains the observed signatures 
of RCW 120, i.e., its ring morphology, coexistence of the two clouds and their 
large velocity separation, and absence of the expanding motion. 


Subject headings: ISM: clouds — Radio lines: ISM 
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Introduction 


Feedback effects from high-mass stars such as intense ultraviolet (UV) radiation, 
strong stellar winds, and supernova explosions at the ends of their lives provide significant 
effects on the evolution of the interstellar medium (ISM) and star formation in galaxies. 
Understanding the formation of high-mass stars is therefore one of the most pressing issues 
in the contemporary astronomy. Despite numerous theoretical and observational studies 


attempting to address this issue in th e last f ew decades, 


involved are still elusive (e.g., 


Porteg ies Zwar t et al 


he detailed physical processes 


20101 ). For instance the gas mass 


accretion rate is one of the known difficulties; high-mass stars must be formed within 


accretion rates of 10 

4 -10 3 M^v r 1 (e.g.. 

Wolfire & Cassinell 

1 | 

1987; 

Tan & McKee 

2003 

Krumholz et al. 

2009; 

Hosokawa & Omukai 

2009; 

Kuiper et al. 

2011 

)• 


Recent attempts to investigate the mechanism of the high-mass star formation in the 
super star clusters (SSCs) Westerlund2 and NGC3603 have provided an important clu e 


to un derstan d the high-mass star formation (jFurukawa et al 


2009 


Ohama et al. 


2010 


Fukui et al. 


2014bh . SSCs have remarkably concentrated stars with total mass of ~ 10 4 M 0 
within a small volume less than lpc 3 , and their origin was not clarified. The authors 
identified a pair of two giant molecular clouds with ~ 10 5 M 0 associated with each SSC, 
where the two clouds has a large velocity separation of 10 - 20 km s -1 . The large velocity 
separations cannot be explained by expanding motion of the gas driven by stellar winds 
or supernovae. They suggested that a collision between the two clouds followed by strong 
shock compression of the molecular gas led to the formation of the massive clusters in a 
short time less than 1 Myr. In addition, a molecular line study in the Galactic Hll region 
Trifid Nebula (M20), which is dominated by a single O star has shown that M20 is a mimic 
of Westerlnd2 and NGC3603, where a collision between two molecular clouds with masses 
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of 10 3 M 0 triggered the formation of the O star flTorii et ah 201l[h Here again, the velocity 
separation of the two clouds (~7kms _1 ) is much larger than the permitted gravitational 
binding velocity of the two clouds. These three cases, Westerlund2, NGC3603, and M20, 
indicate that, although more samples and quantitative studies are necessary, cloud-cloud 
collision has a possibility to be one of the important modes of high-mass star formation for 
a wide mass range, from a single OB star to SSCs. 


ROW 120 is a Galactic Hu region (IRodgers et ahl 196011 located at a distance of 


or CD—38°11636 (F 

'"‘•7 

mssei] 

2003; 

Avedisova & Kondratenko 

1984; 

Georgelin & Georgclin 

1970; 

Zavagno et al 

2007 

). It shows a beautiful ring-like structure in the mid-infrared 


and sub-mm wavelengths which encloses the Hll region. Figure Q] shows a composite color 


emission (green. 

Beniamin et ah 

2003 

). The Spitzer /MIPSGAL 24 mn emission (blue, 

Carev et ah 

2009 

) is distributed inside of the 8 /jm ring, which is attributed to hot dust 


grains illumina ted by t he UV radiation. Far infrared emission at 250 /irn by Herschel /SPIRE 


(red, 


Zavag no et ah 


2010|) also has a ring-like structure, while it is distributed in the slightly 


outer part o f the ; 8 /i rn ring. The age of the exciting star in RCW120 was investigated by 


Martins et ah 


( 20101 ). and the authors discussed that, due to the large uncertainty on the 


effective temperature and the relative low luminosity of the star, any age less than 5 Myr is 
possible. 


RCW 120 is cataloged as the Spitzer bubble S 7 by 


Churchwcll et. ah 


(120061 1. They 


identified about 600 ring-like 8 /im structures within ±60° in the Galactic longitude, 
and th e work was followed by an expanded catalog of 5106 bubbles publish ed by 


Simpson et ah 


Kendrew et ah 


201211 . and many of the bubbles include Hll regions (iDeharveng et a 


1 


20121 ). The pressure-driven expanding Hu region (e.g., 


Weaver et ah 


2010 


1977|; 
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Ho sokawa & Inutsuka 2005 ) is usua 


bubbles (e.g., 


Deharveng et ah 


ly discussed to explain the formation of the Spitzer 


20101) . in which the expanding ionization front preceded by 


the shock front driven by the Hll region accumulates neutral material between the two 
fronts, leading formation of a dense layer. One can observe this dense layer as an expanding 
spherical shell on the sky. The dense layer is fragmented into gravi tationa lly unsta ble 


cond 


ensations, and the condensations are expected to form stars (e.g., 


2001 


Hosokawa fe Inutsuka 


Zayagno et ah 


(120071 ) observed millimeter dust continuum emission in ROW 120 and 


found that ROW 120 is surrounded by dense dust condensations which trace the southern 



careful analysis of infrared point sources using the 2MASS and Spitzer datasets to identify 


the young ste 
observations, 


lar objects ( YS Os) in RCW120 (see Figured]). In addition, with the Herschel 


Zavagno et al. 


( 2010! ) discovered one young massive object identified as Class 


0 in the most massive condensation (condensation lcore) located at the south-west of the 
ring, The stellar mass estimated for the Class 0 object is 8-lOM 0 . Condensation 1 also 


harbors a chain of several Class I anc 

Class I-Class II YSOs, which are aligned parallel 

to the ionization front ( 

Zavagno et al. 

2007; 

Deharveng et al. 

2009 

). These observational 


signatures indicate that RCW 120 is a site of on-going star formation. 


Zayagno et al 


(120071 ) and 


De harven g et al. 


(120091 1 discussed the triggered star 


formation in the ring of the dense neutral medium accumulated by the pressure driven Hll 
region. They suggested that several different star formation mechanisms are simultaneously 
working in the ring, e.g., gravitational instabilities, dynamical instabilities, and interaction 
of the Hll region with the pre-existing molecular condensations. In addition, they discussed 
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that the surrounding ring of RCW120 is porous, and the YSOs located distant from the 
ionization front (e.g., those around the condensation 5 region) could have been formed via 
interaction with the ionized gas leaking through the ring. 

Despite of the detailed studies on the second generation star formation in RCW 120 in 
terms of the triggered star formation, the origin of the exciting O star inside the 8 /im ring 
has never been focused on. In this paper, we present new molecular observations in CO 
<7=1—0 and <7=3-2 toward RCW 120 performed with the NANTEN2, Mopra, and ASTE 
telescopes. Mopra and ASTE provide detailed molecular distributions at high angular 
resolutions, while NANTEN2 is used to describe a large-scale gas distribution, allowing us to 
have the comprehensive understanding of the origin of RCW 120. This paper is organized as 
follows; Section 2 summarizes the observations and Section 3 the results. The discussion is 
given in Section 4, in which two ideas are tested to understand the observational signatures, 
i.e., evolution of the Hu region and cloud-cloud collision. Section 5 gives a summary. 


2. Observations 


The NANTEN2 4-m mm/sub-mm telescope in Atacama, Chile (4,850 m) was used 
to image 1.5° x 1.1° of RCW120 in the 12 CO and 13 CO ,7=1-0 emission by using the 
on-the-fly (OTF) mode during from 2012 May to 2012 December. A 4K cooled SIS 
mixer receiver provided a typical system temperature of ~250 K in DSB, and a digital 
spectrometer provides 16,384 channels at a bandwidth and resolution of 1 GHz and 61kHz, 
which corresponds to 2600kms _1 with a velocity resolution of 0.17kms -1 , respectively, 
at 110 GHz. We smoothed the obtained data to a velocity resolution of 0.75kms _1 and 
angular resolution of 200". The pointing accuracy was confirmed to be better than 15" with 
the daily observations of the Sun and IRC+10216 (R.A., Dec.)(J2000) = (9 h 47 m 575406, 
13°16 ' 43 , . , 56). The absolute intensity calibration was done by daily observations of IRAS 
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16293-2422 in LDN 1689 (R.A., Dec.)(J2000) = (16 h 32 m 23!3, -24° 28' 39"2) and Perseus 
molecular cloud (R.A., Dec.)(J2000) = (3 h 29 m 19!0, 31°24 ' 4970). We compared the 
obtained data with the FCRAO CO data set, where the data were spatially smoothed to 
the same angular resolution with NANTEN2, and estimated the main beam efficiency as 
~0.53. The typical r.rn.s. noise fluctuations in the 12 CO J= 1-0 and 13 CO J=l-0 emission 
are 0.5 K and 0.4 K, respectively. 


The 22-rn ATNF (Australia Telescope National Facility) Mopra mm telescope in 
Australia was used to image RCW120 in 12 CO J=l-0, 13 CO J=l-0 and C 18 0 J= 1-0 
with a high spatial resolution of 33" during 2012 July. We covered 16' x 12' of RCW 120 
in the Galactic coordinate, and the OTF mode was used with a unit field of 4' x 4'. The 
typical system noise temperature, T sys , was 400-600 K in the SSB. The Mopra backend 
system “MOPS” provided 4096 channels across 137.5 MHz in each of the two orthogonal 
polarizations, and the effective velocity resolution was 0.088 km s^ 1 and the velocity coverage 
was 360kms _1 at 115 GHz. The spectra were gridded to a 15" spacing, and smoothed 
to a 45" beam size with a 2D Gaussian function. The pointing accuracy was checked 
every 1 hour to keep within 7" by observations of 86 GHz SiO masers. We used Orion-KL 
(R.A., Dec.) = (—5 h 35 m 14!5, —5°22 , 29 , .'6) for the absolute intensity calib ration and made 


comparisons with the peak temperature of 100 K at Orion-KL given by 


Ladd et al. 


(2005) 


calibrated with an extended beam efficiency introduced in their paper. Consistency of the 
calibrated intensity was then checked by comparing with the NANTEN2 CO J=l-0 dataset. 
We smoothed the channels in velocity to OffikmsW The typical r.rn.s. noise fluctuations in 
the 12 CO, 13 CO and C 18 0 J=l-0 emission are 0.4K, 0.2K and 0.2K, respectively. 


12 COJ = 3 — 2 line observations of RCW120 were performed during 2014 June by using 
the ASTE 10-m telescope situated in Atacama, Chile (4850 m). The SIS mixer receiver 
“CAT345” and the digital spectrometer “MAC” which provided 512 MHz (^ 445 km s" 1 ) 





bandwidth and 0.5 MHz 0.43 kms -1 ) resolution were used. The observations were 
made with the OTF mode at a grid spacing of 7/5 toward the same area with the Mopra 
observations, and the HPBW was 22" at the CO J — 3 — 2 frequency. The pointing accuracy 
was checked every 1-2 hours to be better than 2", and the absolute intensity calibration 
was made by monitoring W28 (R.A., Dec.)(B1950) = (17 h 57 m 26!8, —24°03 / 54'/0), W44 (R.A., 
Dec.) (B1950) = (18 h 50 m 46/l, 1°11'1T.'0), and IRC10216 (R.A., Dec.)(B1950) = (9 h 45 m 14/8, 
3°30 / 40 / . / 0 ) every 1.5 hours. The obtained spectra were compared with the CSO 10 m data 


(Wang et al 


19941 1 and then the main beam efficiency was estimated to be 0.62 ± 0.08. The 


obtained data were smoothed to the angular and velocity resolutions same as the Mopra 
dataset, and the typical r.m.s. noise fluctuations are finally 0.1 K. 


Results 


3.1. Large scale molecular distributions with NANTEN2 


Using the NANTEN2 12 CO dataset, molecular gas distributed toward RCW120 is 
identified mainly in two velocity ranges. Figures [2](a) and (b) show large scale NANTEN2 
12 CO integrated intensity distributions in the two velocity ranges. In the red-shifted velocity 
range at ulsr = — 15-+2kms _1 (Figure Eli), one large molecular cloud which clearly has a 
ring-like component is distributed at ~ 8kms _1 , which shows an excellent agreement with 


the 8 fj ,m ring, indicating physical association with RCW 120. We hereafter re 


cloud as “the red cloud”. By assuming the flat rotation of the Galactic disk flBrand fe Blitz 


er to the large 


19931 ). the distance of the red cloud can be estimated as ~1.3kpc with ulsr of —8kms 1 , 
which corresponds to the previously k nown distance o f RCW 120 estimated with the visual 


extinction towards the central O star flZavagno et al. 
thus estimated to be about 12 pc. 


20071 ). The size of the red cloud is 
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In the blue-shifted velocity range at vlsr — —37—24 kms” 1 shown in Figure HJb) , CO 
distribution is more complicated and localized relative to that in the red-shifted velocity 
range. A CO cloud having a “V-shape” is located at the north of the opening of the 8 /irn 
ring at u L sr of —28 kms 1 . At each side of the V-shape, one CO peak is seen at its lower 
part, and the bottom of the right side of the V-shape looks touching the tips of the 8 /mi 
ring in projection, suggesting a possible interaction with RCW120, although the present 
coarse spatial resolution cannot establish the physical association with RCW 120. We 
hereafter refer to the V-shape cloud as “the blue cloud”. If the blue cloud is associated with 
RCW 120, the size of the V-shape is estimated to be about 11 pc x 7 pc. In the blue-shifted 
velocity range another CO cloud overlapping with the 8 /irn ring is seen around the southern 
part of the ring (7, b S 2 348?20, 0?44) at ~ —33kms -1 (hereafter “the —33kms” 1 cloud”). It 
has weak CO intensity and looks connected with the elongated CO features distributed at 
the southwest of RCW 120 in Figure[2](b). The association of the —33kms” 1 cloud with 
RCW 120 is not clear. 


A comparison between the two velocity ranges is shown in a close-u p in Figure [21(c) , 


on which the dust condensations and YSOs identified by 


Deharveng et al. 


< 20091) are 


superimposed. The CO peak on the right side of the V-shape looks coinciding with the 
local CO peak in the red cloud at (/, b) ±2 (348?33, 0?64), while the CO peak on the other 
side seems to have complementary distribution with the rim of the red cloud. YSOs are 
extensively distributed throughout the red cloud. Class I and Class I-Class II YSOs (red and 
white circles, respectively) are mainly distributed around the molecular ring component in 
the red cloud, in which many dust condensations are also concentrated, while Class II YSOs 
(white triangles) are mainly seen toward the outer part of the red cloud. Condensations 6 
and 9 are distributed toward the coinciding peaks of the blue cloud and the red cloud at 
the north of the 8 /irn ring, and several YSOs are also seen there. 
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Position-velocity diagrams for a large area of RCW120 are presented in Figure [3] 
with CO (contours') and Hi (i mage), w here the SGPS(Southern Galactic Plane Survey) 


Hi dataset is used (IMcClure-Griffiths et al. 


20051) As a supplement, CO velocity channel 


maps including the Mopra datasets are shown in Figures IT5l IT9l in Appendix. The CO 
emission in the red cloud and the blue cloud have typical velocity widths of ~10kms _1 and 
~5kms _1 , respectively, and show almost uniform velocity distributions well beyond the 
extent of RCW 120. In Figure[3]the Hi emission in the red cloud is seen as an intensity gap 
in the bright Hi emission extended over ±15kms _1 , while that in the blue cloud is seen at 
around —26kms _1 , which is slightly shifted from that in CO around — 28kms~ 1 . 

As depicted by the red arrow in Figure[3l the blue cloud has a velocity feature at 
(/, b) ^ (348?25, 0?65) extended to ~ — 22kms -1 toward the positive velocity side. The 
spatial distribution of the high velocity feature is shown in Figure[21(c) with the thick black 
contour. It has a size of ~ 0?4, which corresponds to ~ lpc at the distance of RCW 120, 
and is located just next to the coinciding peaks of the blue cloud and the red cloud. The 
corresponding feature is also seen in the Hi emission (Figure[3|). The Hi feature has a much 
larger velocity extent and appears to be connecting the blue cloud with the red cloud, 
suggesting a physical interaction between the two clouds and RCW 120. In addition, Hi 
has other several high velocity features throughout the blue cloud as shown by the blue 
arrows with the dashed lines in Figure[3] These Hi features are not seen at the negative 
velocity range, and provides another support for the physical connection between the red 
cloud and the blue cloud even outside RCW 120. The typical atomic column density of the 
high velocity Hi features is ~ 4 x 10 20 cm -2 by assuming optically thin 21-cm emission, and 
is 
in 


~ 8 x 10 20 cm 2 bv considi 

Fukui et al. 

(2014a, 

2015a) 


features is estimated to be ~ 100 1 
adopting the opacity correction by 


Tm for the optically thin ca se, and to be ~ 200 M 0 by 


Fukui et al 


(2014a, 


2015a) 
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3.2. Detailed molecular distributions with Mopra 


Detailed CO distributions toward RCW120 using the Mopra datasets are shown 
at 0.2 pc spatial resolution. Figure ® shows the integrated intensity distributions of the 
12 CO, 13 CO and C 18 0 emission in the red-shifted and blue-shifted velocity ranges. In 
FiguresH)(a)-(c), the red cloud clearly shows a ring-like structure in all the three CO lines, 
where the molecular ring consists of small clumps with sizes of 0.3-0.5 pc. Detection of 
C 18 0 suggests the existence of the dense gas in th e molecular ring, and the dense molecu lar 


ring corresponds well with the dust distribution (jZayagnp et al 


2007t 


De harven g et a l. 


20091 1. Typical H 2 column density 1V(H 2 ) of the ring component is estimated to be 


20 


~3x 10 22 cm -2 with the 12 CO i ntegra ted intensities by assuming an X-factor of 2 x 10 
(Kkms -1 ) -1 cm -2 ((Strong et al.l 1198811 . where the X-factor is an empirical conversion 
factor from the 12 CO integrated intensity to 1V(H 2 ). On the other hand, 7V(H 2 ) toward 
the O star located inside the ring is estimated to be 1.4 x 10 22 cm -2 . This estimate can 
be compared with the visual exti nction A v derived as 4.65 mag toward the same direction 


(jAvedisova & Kon drate nko 


1984| ). which corresponds to the total ISM column density 
A r (Hi+H 2 ) of 9 x 10 21 cm -2 , almost comparable to the above column densities, by u sing the 


relati ons -/VfHl+Fhd = 5.8 x 10 21 E(B — V) cm 2 mag 1 and A, = 3.1 E(B — V) ( Bohlin et al. 


19781 ). It is thus suggested that the exciting O star is on the far side of the red cloud along 


the line-of-sight. 


CO distributions in the blue-shifted velocity range obtained with Mopra are shown 
in Figures[U(d)-(f). The blue cloud is detected in C 18 0 with a more than 3 a level, while 
the — SSkrns^ 1 cloud is seen only in 12 CO at the southern part of the 8 /irn ring. It is 
notable that the blue cloud has a molecular protrusion just next to the western tip of the 
8 /mi ring, and it looks elongated inside the 8 /irn ring through the opening at the north. 
Typical 1V(H 2 ) of the «»33kms -1 cloud and the blue cloud are estimated with 12 CO to be 
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~ 4 x 10 21 cm 2 and ~ 8 x 10 21 cm 2 , respectively, significantly smaller than that of the 
ring component of the red cloud. 


In Figure[5] comparisons between CO and other wavelengths are presented. FiguresEfa)- 
(c) show the 13 CO integrated intensity distributions of the red cloud (black co ntours and 


colors) superimposed on (a) the Spitzer 8 pm, (b) the SuperCO SMOS Ha flParker et al. 


200511 . and (c) the 843 MHz radio continuum flBock et ah 


19991 1 grayscale images. In the 


SUMSS image in Figure[5](c), the thick contours are plotted at the 5 a level. As seen in 
Figure [51(a) , the 8 pm ring apparently traces the inner rim of the red cloud, supporting 
physical association between the red cloud and RCW120. No apparent molecular 
coimterpart is seen toward the western tip of the 8 pm ring, which is also confirmed with 
the Herschel 250 pm emission in Figure]!] In the supplemental velocity channel maps in 
FigureflSl a possible counterpart for the western tip is seen in 12 CO at around —lOkms -1 , 
although it is weak in 13 CO and not detected in C 18 0. Distribution of the ionized gas 
traced with Ha in (b) and the radio continuum emission in (c) is enclosed by the 8 pm ring. 
The filamentary dark lanes running across the Ha image are due to extinction and has CO 
counterparts in the red cloud, indicating that these molecular counterparts are located on 
the near side of the Hu region. 


In Figures E(d)-(f) 12 CO distribution of the blue cloud is used to make comparisons 
as done for the red cloud in Figures[5](a)-(c). As also shown in Figure|U(d), a molecular 
protrusion is elongated from the blue cloud toward inside the 8 pm ring through its opening, 
while the —33kms _1 cloud having very weak intensities shows no clear correlation with 
the infrared emission. Although the spatial correlation between the blue cloud and the 
8 pm emission is not apparently identified, the ionized gas clearly shows complementary 
distribution with the blue cloud, especially with the molecular protrusion, as seen in 
FiguresEte) and (f), where the lowest contours are plotted at the 5 a level in the radio 
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continuum image. Since the complementary distribution is seen not only in Ha but also in 
the radio continuum, it is not due to the extinction but is attributed to exclusion of the 
ionized gas by the blue cloud. This result suggests that the blue cloud as well as the red 
cloud is possibly associated with RCW120. 


3.3. CO J— 3 — 2 / J— 1—0 intensity ratios and LVG analysis 


In order to investigate the association of the molecular gas with RCW 120 into 
more detail, distributions of the 12 CO J=3-2/J=l-0 intensity ratio (hereafter i? 3 _ 2 /i_o) 
are presented with the Mopra and ASTE datasets. Intensity ratios between different J 


(see e.g., 


Torii et ah 

2011; 

Fukui et ah 

2014b1 


2014bl ). i? 3 _ 2 /i—o distributions of the red cloud and 


the blue cloud are shown in the upper two panels in Figured while in the lower two panels 
the region with high i? 3 _ 2 /i_o > 0.8 is presented in contours superimposed on a composite 
image of the Herschel 250 /an (red) and the Spitzer 8/mi (green). R 3 _ 2 / 1 _ 0 typically ranges 


from 0.6 to > 1.2 in both of the two clouds. Considering t 


rat typica. 


clouds without any star formation have R 3 - 2 / 1-0 of ~ 0.4 ( Oka et ah 
present two clouds is significantly high. 


i? 3 - 2 /i-o in the local 


20071), -R 3 - 2 / 1-0 hi the 


I 11 the red cloud, high R 3 _ 2 /i-o of > 0.8 is mainly seen where the dense molecular 
ring is located, which is also traced with the 250 pm emission (FigurelHt). The highest 
-R 3 - 2 / 1-0 of > 1.2 is seen around the southern part of the ring, in which condensations 
1-3 are distributed. Another high i? 3 _ 2 / 1 „ 0 region is seen inside the 8/mi ring at 
M) (348?30, 0?50-0?55), whose western rim is traced by the 8/mi and 250 p emission, 
suggesting that the cloud is illuminated by the O star. The western tip of the 8 /mi ring 
in contrast shows relatively low R 3 _ 2 /i_ 0 compared with the other part of the 8 /mi ring. 
On the other hand, i? 3 _ 2 /i_o in the blue cloud is high at its southern part. The highest 
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i? 3 - 2 /i-o is seen at an extension of the eastern tip of the 8 /im ring, while the western 
tip also shows high i? 3 _ 2 /i_o- As seen in Figure^d), these high i? 3 _ 2 /i_o regions show 
excellent coincidence with the 250 /im emission, particularly for the very tip of the western 
side of the 8/an ring. Another high R 3 _ 2 / 3 _ 0 is distributed toward condensation 5a at the 
south-eastern rim of the blue cloud, suggesting possible association between the blue cloud 
and the condensation, although the red cloud also shows high R 3 _ 2 /i_ 0 toward the same 
direction as shown in Figure^a). 


In order to investigate the origin of the highly excited gas in both the red cloud 


and the blue clouc 


Goldreich fe Kwan 


, we here perform the large velocity gradient (LVG) analysis (e.g., 


1 9741 ) to estimate the kinetic temperature T ^ and number density n(H 2 ) 


of the molecular gas, using the intensity ratios 13 CO J=1-0/ 12 CO J=l-0 (hereafter R\ 3 /i 2 ) 
as well as R 3 - 2 / i_o- It should be noted that the assumption of the uniform velocity gradient 
is not always valid in the molecular gas associated with Hu regions, but radiative transfer 


calculations considering a microturbulent cloud in 
significant difference from the LVG analysis(e.g., 


eracting with an H l l region show s no 


Leung fc Liszt 

1976; 

White 

1977) 


therefore adopt the LVG approximation in the present study. 


Six target regions (regions A-F) used in the present analysis are shown in Figures|6](a) 
and (b). Since 12 CO is strongly affected by self-absorption toward the dense part of the 
ring, particularly around condensations 1 — 3, these regions are not adopted. In the red 
cloud regions A and B are distributed toward condensations 4 and 5c in the dense ring, 
respectively, while region C is chosen outside the dense ring to make comparisons with 
regions A and B. In the blue cloud all the three regions are chosen toward the high R 3 _ 2 n_ 0 
regions around the opening of the 8 /irn ring. 


Details of the setup of t he ana lysis is_as follo ws: we adop 
[ 12 CO]/[H 2 ] = 10- 4 (e.g., 


Frerking ct al. 


1982 


Leung et al 


an a bundance ratio X(CO) 


19841) and [ 12 CO]/[ 13 CO] 
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= 77 (WUson & Rood 


1994]), and assume uniform velocity gradient dv/dr in each of the 


red cloud and the blue cloud as 6kms _ 1 /0.3pc ~ 16.7kms -1 pc -1 and 3kms _ 1 /0.3pc 
~ 6.7kms _1 pc -1 , respectively. Errors in i? 3 _ 2 / 1_ 0 and i ? 13 / 12 are calculated with la 
baseline fluctuations of the spectra and 10 % relative calibration error between two CO 
transitions. The 10% relative calibration error is applied only to i? 3 _ 2 /i_ 0 . Since 12 CO 
J = l -0 and 13 CO J = l-0 were obtained simultaneously with the same receiver and the 
same backend system at Mopra, the relative uncertainties between these two transitions are 
expected to be negligibly small. 


Figurc[7] shows the results of the LVG analysis, in which T ^ and n(H 2 ) are determined 
where the i? 3 _ 2 /i_o area (blue) and the -R 13/12 area (red) overlap. Regions A and B at the 
ring component show high temperature of > 40 K at high density > 10 5 cm -3 , while region 
C outside the ring indicate relatively low temperature of ~ 20 K at 10 4 cm -3 . The blue 
cloud also shows high temperature similar to those in regions A and B, but seems to have 
lower n(H 2 ) of > (1-4) x 10 4 cm -3 . 


Possible heating mechanism includes the following two; one is heating by embedded 
stars in the molecular gas. The high 21 gas with counterparts in YSOs or dust condensations 
may be understood with this heating, which may be the case in regions A, B, and D. 

The other is heating by the feedbacks from the exciting O star, i.e., UV radiation and/or 
interaction with the Hll region. This heating probably works in the blue cloud, as shown 
with a negative i? 3 _ 2 / 1-0 gradient from the south to the north seen in the blue cloud 
(FigureEb), suggesting that it depends on the distance from the Hll region in RCW120. 
The high R 3 _ 2 /i_o region which seems to be in contact with the Hll region (see Figure[5]) 
may be interpreted as due to the interaction with the Hll region. On the other hand, the 
region around the 8 pm PAH ring may be understood with the radiative heating by the O 
star. Region E may be a case of the interaction with the Hll region, while region F keeps 
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both possibilities, and the high i? 3 _ 2 /i_o region at the western tip of the 8 /ini ring in the 
bine cloud may be explained with the radiative heating by the O star. 

As a summary, the physical association of the red cloud and the blue cloud in RCW120 
is strongly indicated by the following three results; 

1. High .R 3 _ 2 /i— o and high kinetic temperature (>30-40 K) of the southern part of the 
blue cloud (Figures^)] and [TJ, which coincides with the 250 gm image, indicate the 
heating by the feedbacks of the exciting O star in RCW 120. 

2. The complementary distribution between the blue cloud and the ionized gas tracers 
(i.e., Ha and 845 MHz emission) supports physical interaction between them 
(Figure[5]). 

3. The intermediate velocity features which connect the red cloud and the blue cloud are 
seen especially in the Hi emission (Figure[3]). 


3.4. Mass of the molecular gas in RCW 120 


By adopting the same distance of 1.3 kpc both for the red cloud and the blue cloud, 
their molecular masses can be estimated by using the integ rated intensity of t he 12 CO 


emission with an X-factor of 2 x 10 20 (K km s 3 ) 1 cm 2 (Strong et al.l 1 9881 ). The 


molecular mass of the entire red cloud shown in Figure [21(a) is calculated to be 5.1 x 10 4 M 0 
with the NANTEN2 12 CO dataset, and that of the ring component is to be 4.2 x 10 3 M 0 
with the Mopra dataset, where the target region is shown by thick black contours plotted 
at 100 K kms -1 in Figure IU(a). On the other hand, the molecular mass of the blue cloud in 
Figure [2](b) is calculated to be 7.5 x 1O 3 M 0 . 
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3.5. Physical condition of the ionized gas in RCW 120 


To estimate the mass of the ionized gas by measuring the electron density n e , the 
average flux densities of the 843 MHz emission (FiguresEfc and f) is here measured for the 
southern half and the northern half of the Hll region as separated by the dashed line in 

ure T 0 of 8000 K an d thepath length 


Wood fe Churchwell (119891 ) enable us 


Figure[5](f). By assuming uniform electron tempera: 

L of the Hu region of 3.5 pc, the equations 2-4 in 
to calculate n e to be 110 cm -3 and 60 cm -3 for the southern half and northern half of the 
Hu region, respectively. If T e is varying from 5000 K to 10000 K, the results change only by 
±10%. The total mass of the ionized gas including ions is thus estimated to be ~ 140 M 0 
by assuming a spherical shape of the Hll region. 


The derived n e is used to estimate the pressure of the Hu region as P e = 2 n e x T e = 
1.8 x 10 6 Kcm -3 for the southern region and 10 6 Kcm -3 for the northern region. These 
figures can be compared with the turbulent pressure of the surrounding cold neutral 
medium P n . The molecular number densities n (IF) is given in the LVG analysis (Figure!?]) 
as 10 5 cm -3 for the red cloud and 10 4 cm -3 for the blue cloud. On the other hand, by 
using the linewidths of the 13 CO profiles, effective temperature of the turbulent gas can be 
estimated to be 550 K and 140 K for the red cloud and the blue cloud, respectively. P n in 
the red cloud is finally given as 5.5 x 10' Kcm -3 , one order of magnitude larger than P Cl 
while 1.4 x 10 6 Kcm -3 in the blue cloud, which is nearly equal to P c . 


These estimates indicate that the Hll region cannot drive the bulk motion of the red 


cloud, while th e effect may be margina 


as discussed in 


Walch ct al 


(20U 


for the blue cloud. It is, however, noted that, 


20121 ), the ionized gas can stream through less dense 


gas if the cloud has an inhomogeneous density distributions, and this may be the case in 
RCW 120. As shown in Figures[5l and f, the distribution of the ionized gas coincides with 
a gap between the two dense clumps in the blue cloud, implying erosion by the Hll region. 
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On the other hand, the dark 8 /mi feature seen in north-east of the 8 //in ring (Figure[5]) 
has many filamentar y su bstructures which are ra dially elongated. For these structures, 


Zavag no et al. 


(120071 ) and 


De harven g et al. 


(120091 ) discussed leaks of the ionized gas through 


the dense ring component. The leakage of the ionized gas may l ead to un deres timate of 


the total mass ionized by the O star. However, as discussed by 


Anderson et al 


(120141 ). the 


mass fraction of the escaping ionized gas may be minor, about 25 % of that observed. 


4. Discussion 

Our new CO observations with NANTEN2 and Mopra have shown that the two CO 
clouds, the red cloud and the blue cloud, at velocities of ~ —28kms _1 and ~ —8kms _1 are 
both associated with RCW120. The red cloud with a mass of 5 x 10 4 M 0 clearly traces 
the 8 //m ring, while the blue cloud with a mass ~ 8 x 10 3 M 0 located just north of the 
opening of the 8 /mi ring shows complementary distribution with the Hll region. Both of 
the two clouds have high kinetic temperature > 20-40K, supporting the coexistence of the 
two clouds at RCW 120. 

To understand the coexistence of the two clouds, interpreting the large velocity 
separation of ~20kms _1 may be the key. By tentatively assuming an viewing angle of 
the clouds’ relative motion to the line-of-sight as 45°, the relative velocity between the 
two clouds is estimated to be ~30kms _1 , which requires the total mass of ~ 2 x 10 6 M 0 
to gravitationally bind the two clouds within 12 pc of the exciting star. This figure is 
significantly larger than the total molecular mass in RCW 120 derived as ~ 6 x 1O 4 M 0 , 
indicating that the coexistence of the two clouds cannot be understood as a gravitationally 
bound system. Even if we adopt the projected velocity separation 20kms _1 as the real 
relative velocity, the cloud motion is not gravitationally bound. 
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4.1. Evolution of the Hu region as the origin of RCW 120 

A possible explanation of the velocity separation between the two clouds is the feedback 
effects of the exciting star, such as the stellar wind, UV radiation, and pressure driven 
Hu region. Supernova explosions are also an alternative candidate of the driving source, 
whereas, in RCW 120 no supernova remnants have been identified to date. In the following, 
in order to explain the morphology of RCW 120 and coexistence of the two clouds, we 
examine several ideas in terms of the evolutional scenario of the Hll region. 


4-1.1. Cases of uniform ambient medium 


We first investigate the evolution of the Hll region as the origin of the two clouds in 
RCW 120, by testing whether it can explain the ring-like morphology of RCW 120, velocity 
distribution of the molecular gas, and the morphology of the two clouds. A basic idea 


to explain the ring shape is k nown a s “Co 


first proposed by 


lect & Collapse (hereafter C&C)”, which was 


E lme green & Lada (1977). An expanding Hll region sweeps up uniform 


ambient material between the ionization front and the preceding shock 


dense molecular layer (e.g., 
layer fragments via gravita 


Hosokawa fe Inutsuka 


2005 


2006 


Dale et al. 


'ront to form a 


20071 ). The dense 


a new generation of stars flWhitworth et al. 


ional instability, and th e i ndividual fragments may finally fo rm 


1994a 

b; 

Elmegreen 

1994; 

Dale et al. 

2009) 


This model is usually d iscu ssed as the origin of the Spitzer bubbles including RCW 120 


(e-g-, 


Deharveng et al. 


20101 ). 


The molecular masses accumulated within the dense layer by C&C, M ring , is estimated 
with a equation 4/37Trf ing po, where r ring is the radius of the ring, and po is the initial density 
of the ambient medium which corresponds to number density no- By assuming M r - mz of 
3100 M q derived in Section 3.5, n 0 is calculated to be ~3000cm~ 3 , which is consistent with 
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the estimate of 


Zavagno et ah 


(120071 ). 


Hosokawa & In uts uka (2005 


20061) provide a basic ID 


hydro dynamical numerical simulations of C&C, in which the radiative transfer o f UV and 
far-UV radiation is taken into account. The case of the model S19 in 


Hosokawa fe Inutsuka 


(120061) . in which the mass of the exciting star is 19.0 M 0 , provides nearly the same situation 


of RCW120, and by using their equation (31) and (36), the formation timescale of RCW120 
with a rin g radius o f 1.7 pc can be estimated to be ~0.4Myr (see Figures 14 and 15 in 


Zava gno et al. 


20071) 


A schematic of the characteristic velocity distribution predicted by C&C is shown in 
Figure [HI If the dense layer surrounding the Hll region expands at an expanding velocity 
Wexp because of high pressure driven by the ionized gas, the expanding spherical layer is 
expected to have velocity distributions like those shown in Figures [8](b) - (c). If velocity 
width, An, is small enough, the near and far sides of the layer can be separately observed 
toward the exciting star with a velocity separation of 2n exp (Figure[8}j). On the other hand, 
if An is too large compared with n exp , the near and far sides cannot be distinguished, and 
only one single component with an elliptical shape may be observed, where the direction of 
the O star has the largest velocity width (Figured). 


Figured shows the longitude-velocity map of the Mopra 12 CO data at a latitude range 

le NAN T EN2 datas et in FigureEJ n exp at 


0?42-0?58, which is complementary to that of t 
the time 0.4 Myr is estimated to be ~2kms -1 flZavagno et al 


20071 ). We here test two 


cases; one is the case that both of the red cloud and the blue cloud were formed by C&C 
(case I), and the other is that only the red cloud which shows a ring-like shape was formed 
through C&C, and the blue cloud is associated with RCW 120 by some other reason (case 

II). 


For the case I, as clearly seen in Figures[3] and [9](a), both of the red cloud and the 
blue cloud have uniform velocity distributions well beyond the extent of RCW 120, and 
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furthermore, the velocity separation of ~20kms” 1 is much larger than the expected velocity 
separation of 4 kms” 1 in the C&C scenario. To discuss the case II, 13 CO and C 18 0 velocity 
distributions of the red cloud is shown in Figure [21(b) ■ Our datasets indicate only one 
velocity component toward the exciting star, and the velocity width toward the exciting 
star is smaller than those in the other parts, showing no signature of the expanding motion. 


In order to show the velocity distributions into more detail, intensity-weighted 
velocity map and velocity dispersion map are shown by the first-moment map and the 
second-moment map, respectively, in FiguresdOl Although velocities are slightly shifted 
from the ring component (~ —8 kms” 1 ) to the ring interior (~ —8.Skins” 1 ) as seen in 
Figure fTOf ah the velocity dispersion shown in Figure fTOl b) is increased toward the ring 
component (~ 2-2.5kms” 1 ) relative to that toward the ring interior and the exciting O 
star 1.5kms” 1 ), which is an opposite trend to what is expected in the C&C model. 


If the Hll region is in the pressure equilibrium with the surrounding neutral 
medium, which is possible as discussed in Section 3.5, the expanding motion may not be 


observable, and the o 


bserva 


Hosokawa fe Inutsuka 


ional signatures may become consistent with the C&C model. 


(120061 1 provide the relationship between the expanding velocity and 


the time in their equation 44. To reduce n exp to be less than our detection limit of 0.5kms 


-l 


at the gi ven rip of_3000cm 3 in 5Myr, which is the safe upper limit of the age of RCW120 


([Mart i ns et ah 


2010 |), t 


re UV flux of the exci 


(see the equation 40 of 


Hosokawa fe Inutsuka 


ingjstar is required to be < 10 45 ' 2 photons s 


,-i 


20061). It is two orders of magnitude smaller 


than that of the exciting star in RCW 120 (e.g., 


Zavagno et al 


2007]), and it is therefore 


difficult to explain the observations with the pressure equilibrium assumption. 
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4-1-2. Cases of non-uniform ambient medium 


Recent theoretical studies 
ambient medium. Walch et al 


l ave investigated evolution of Hll regions in non-uniform 


(20U 


2012 ! ) presented numerical calculations of the evolution 


of an_HlI region in a spherical cloud having fractal density structures. In particular, in 


Walch et al. 


(1201111 . the authors investigated the formation of RCW120 and successfully 


reproduced the morphology, size, and mass of the neutral medium in RCW 120. In their 
calculations, the ionized gas streams out through less dense regions rather than dense 
regions, and the pressurized low-density region compresses the remaining gas, leading to the 
formation of dense clumps in the cloud. The resultant density distribution strongly reflects 
the initial condition. In their study, the shell-like structure is still forme d by expans ion, and 


the cold shell material is continuously accelerated up to 4kms 1 (IWalch et al. 


2012). The 


discussion for case I and case II with the uniform medium case holds valid; no expanding 
motion having a velocity separation of 8kms -1 is seen in RCW 120. 


Ochsendorf et al. 


(2 0141 1 calculated evolution of an Hll region by assuming a Bonnor- 


Ebert sphere, in which the exciting star has an offset from the center of the sphere. In 
their calculations, first C&C works to form a dense layer, but the layer is eventually broken 
by the over-pressured ionized gas at the thinnest point of the layer. The ionized gas then 
streams into the surrounding ISM, forming the observed morphology of RCW 120. This 
scenario can explain the offset position of the exciting star from the center of the 8 /mi ring, 
as well as the morphology of RCW 120. The formation timescale they estimated is ~3 Myrs, 
also consistent with the observational limit (< 5Myrs). To make a comparison with their 
model, azimuthal distributions of the NANTEN2 CO maps are presented in FigurefTTl 
where two inner and outer ellipses comprise the target ring as shown in Figures fTTT a) and 
(b). The inner ellipse is determined by fit to trace the inner rim of the 8 /mi ring, and the 
outer one is plotted to have a 1.5 times larger size to fully cover the strong CO emission. 
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The position 
presented by 


angle of t he e llipse is 47° in a counterclockwise direction. In the scenario 


Ochsendorf et al. 


(]2014| ), CO emission should be strongest at the point where 


the distance from the exciting star to the ring is smallest, which here corresponds to Az = 
0°, and should be weakest around the opening of the 8/im ring (Az ~ 180°). The observed 
CO emission is, however, strongest at Az ~ 100°, and weakest at Az ~ —150°, which is 
inconsistent with their scenario. 


Here is a brief summary of the above discussion. In the previous studies, interpretations 
of RCW120 based on the evolution of an Hii region are basically successful in reproduction 
of the morphology and mass of RCW 120. However, the velocity distribution revealed in 
the present observations shows no expanding motion contrary to what the model predicts. 


4.2. An alternative: Cloud-Cloud Collision 

We here postulate an alternative idea that an accidental collision between the red 
cloud and the blue cloud triggered the formation of the O star and the surrounding ring 
component in RCW 120. The idea provides an explanation for the coexistence of the two 
clouds having a large velocity separation as well as the other issues laid on this object. 


4-2.1. Previous studies on cloud-cloud collision 


There has been discussion that collisions between the interstellar molecular clouds 
may trigger the star formation (e.g., 


Loren 


1976 


Higuchi et al. 


2010 


Duarte-Cabral et al. 


2 0111) . Many of the previous studies reported discoveries of pairs of molecular clouds having 
velocity separations of a few kms -1 . Since the small velocity separation can be understood 
not only as Cloud-Cloud Collision (CCC) but also as the gravitationally bound system, 
these works do not provide compelling evidence for the CCC. Another possible star forming 
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site induced by CCC was found in the Galactic center (IHasegawa et al.l Il994h . This is a 
case having large relative velocity in the collision, ~30kms -1 , with triggered formation 
of O stars and the gravitational binding is unlikely, making the CCC as a plausible 
interpretation. This however does not provide a case of CCC-triggered star formation 
common throughout the Galactic disk, since the physical conditions in the Galactic center 
are significantly different, including highly turbulent state, from that in the Galactic disk. It 
is only recent that robust evidence for the CCC and its role on the high-mass star formation 
has been presented observationally in the Galactic disk, even the number of the case is still 
limited. In the observational studies on Westerlund2, NGC3603, and M20, associations of 


(Furukawa et al. 

2009; 

Ohama et al. 

2010 ; 

Fukui et al. 

A 

T—1 

o 

CN 

Torii et al. 

2011 ) 


velocity separations are too large to be explained by the gravitational motion of a bound 
orbit, and an accidental collision is strongly suggested. Recent finding of another high-mass 
star formation triggered by cloud-cloud collision in the N159W with ALMA len ds su pport 


for the important role of CCC even in the Large Magellanic Cloud (jFukui et al. 


2015b ). 


by 


A pioneering numer ical s tudy of CCC was carri ed out by 


Anathpindika (120101) and 


Takahira et al. 


H abe fe Ohtal ( 1992 ). followed 


(2 0141) . These calculations indicate that 
CCC works as a trigger to form massive molecular clumps whose masses are larger than 


their Jeans mass. In addition, 


Inoue & Fukui (120131 ) presented 3D magnetohydro dynamic 


calculations of CCC and discussed that the magnetic fields can help to form massive 
molecular filaments and clumps by driving gas flow to convex points of the deformed shock 
wave, leading large mass accretion rates of 10~ 4 -1 0~ 3 M ,m yr~L S u ch co llisio ns are expecte d 


to be frequent in gas-rich galaxies as discussed by 


mm- 


Tasker fc Tan] ( 20091) and 


Dobbs et al. 
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4-2.2. Basic scenario of the Habe & Ohta model 


In order to understand the present observations in RCW120, we first discuss the basic 
idea of the CCC scenario presented by lHabe fe O hta (119921 1. in which collisions between 


two molecular clouds having different sizes were investigated. Because the authors did not 
include star formation and subsequent UV radiation, in the schematic shown in Figure 
[i~2l some relevant events are qualitatively added. First, a small cloud and a large cloud 
are approaching each other (stage 0). Once a collision occurs, a dense layer is formed by 
compression at the interface of the two clouds, creating a cavity in the large cloud (stage 
1). The diameter of the cavity is nearly equal to that of the small cloud. The compressed 


layer is highly 


(lHaworth. et al. 


urbu 


ent, which can be observed as an intermediate velocity structure 


2015 ). and the small cloud streams into the compressed layer during 


the collision. Along with the cavity creation, massive clump formation and subsequent 
high-mass star(s) formation occur in the compressed layer, followed by the strong UV 
radiation from the high-mass star(s) which ionizes the surrounding neutral material (stage 
2 ). At the final stage (stage 3), the ionized gas fills up the cavity and erodes the inner 
surface of the cavity which includes the compressed layer, enhancing the PAH emission 
observable at 8 /irn. 


Stage 3 may have variations depending on the velocity of the compressed layer. For 
example, if deceleration is fast, cavity creation almost stops and much molecular gas 
remains forward of the moving direction (stage 3A). On the other hand, if the velocity 
keeps fast, the dense layer completely penetrates the large cloud, forming a tunnel instead 
of a cavity (stage 3B). In either case, if the ionized gas completely escapes, only the cavity 
or the tunnel is observable. The depth of the cavity and the location of the O star in the 
cavity are determined by density of the clouds and physical lengths of the small cloud along 
the colliding direction. A detailed analysis of the cavity creation and its effects on the star 
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formation will be presented in a forthcoming paper (Torii et al. 2015b, in prep.). 

The resultant cavity/tunnel should be observed as a projection on the sky. This is 
just like looking at a “wine glass” (stage 3A) or a “tube” (stage 3B) from various angles. 
If observing from the bottom or the top of the glass/tube, one can see a ring shape, and if 
from the sides, one may see a U-shape or two separated components as shown in Figured^] 
as 3A and 3B, respectively. 


It should be noted that, in the Habe & Ohta model it is assumed that two molecular 


clouds with different sizes collide. On t 


have almost the same sizes (ITorii et al. 


re ot 


ler hand, in M20 the observed colliding clouds 


20111 ). and in such a case creation of the cavity is 


not expected. It is likely to be more frequent that collisions take place between clouds with 
different sizes, and the case in M20 is thus relatively rare. 


4-2.3. Cloud-cloud collision in RCW120 

We here discuss an application of the basic CCC scenario to the observed RCW 120, 
where the blue cloud corresponds to the small cloud in Figured^] and the red cloud to the 
large cloud. Figure 1T31 shows an interpretation of the CCC in RCW 120, where the X- and 
Y-axes are respectively taken from the —90°-+90° line and the 0°-180° line defined in 
Figure fTIT a). and the Z-axis is along the line-of-sight. The origin of the new coordinate 
system is taken at the exciting O star. As shown in Figure[2](b), the blue cloud has a 
V-shape structure, each side of which contains one single CO peak at its southern part. We 
assume that there was the third peak with a similar size at the south of the present blue 
cloud and it collided with the red cloud, forming the cavity observed as a ring shape in 
RCW 120. The viewing angle of the collision is assumed to be 45° toward the north from 
the line-of-sight (= Z-axis). The colliding velocity is then calculated to be ~30kms -1 , and 
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the depth of the cavity is estimated to be roughly 5 pc with the observed size of the 8 pm 
ring, ~3.5 pc. 0 star will be formed in the compressed interface 
is accreted in the high-mass dense cores which are simulated by 


layer when su: 

Icient. 

Inoue & Fukui 

(2013) 


interface layer between the two clouds will move at a velocity somewhat less than 30 km s 1 
d ue to th e dec eleration depending on the density ratio between the two clouds (e.g., 


Takahira et ah 


2014( 1. If we tentatively assume that the decelerated velocity is 15kms 1 , 


the O star will move at the same velocity in the same direction with the small cloud. The 
observed position of the O star with an offset from the center of the ring to the southwest 
on the symmetry axis of the cavity is then a natural consequence of the triggering under 
the cloud motion to the southwest. 


After forming the O star in the shocked layer, the UV feedback immediately ionized the 
surrounding neutral material and prohibited formation of other stars via ionization of the 
neutral gas, followed by the rapid evolution of the Hll region. The Hll region finally filled 
out the cavity as observed, indicating that the CCC in RCW120 is currently in stage 3 
(Figured!]). Because much gas still remains toward inside the 8/mi ring, the case 3A is here 
preferable rather than the case 3B. In Figured!] the CO integrated intensity distributions 
of the red cloud along the X- and Y-axes are presented. The plots show significant CO 
detection inside the 8/irn ring both in 12 CO and 13 CO, although it is strongly reduced to be 
40-80% in 12 CO and to be 25-50% in 13 CO from the ring component. In addition, as 
discussed in section 3.2, the O star is likely on the far side of the red cloud, which offers 
another support for the case 3A. 

Differently from the basic model in Figured!! in which at stage 3 the ionized gas 
seems to easily stream out through the opening of the cavity, the opening of the cavity in 
RCW 120 looks being blocked by the bottom of the blue cloud (Figured!]). The column 
density of the blue cloud is by one order of magnitude smaller than the red cloud, making 








28 


it easier for the ionized gas to stream out through it. The Ha and radio continuum images 
in Figure^ however, do not s how a s i g n of such streams. Although diffuse Ha emission is 


seen outside the 8/iin ring ( Deharveng et a h 


20091) . it is very weak and is extended only by 


~l-1.5pc from the 8 /.im ring even around the opening, not showing evidence for the strong 
ionized stream. It is thus s u ggest ed that we are observing the beginning of such event, as 


discussed in 


Zavagno et al. 


(2007). The beautiful symmetric ring struc 


the discussion; as presented in the numerical calculations in 


ring struct 

ure alf 

30 supports 

Dale et al. 

(2012 

2013 

), the 


neutral medium in a mature Hll region may show an asymmetric shape depending on its 
initial density distribution, suggesting that RCW120 is in its early evolutional phase. 


The plots in Figure[[4] are also used to estimate the mass of the molecular gas 
accumulated by the collision to form the cavity. By measuring the intensity gap between 
inside and outside the molecular ring in 12 CO, the accumulated molecular mass is roughly 
estimated to be 1000- 3000 M & . On the other hand, the accumulated layer could have 
included the mass of the partner cloud of the collision. By assuming the H 2 column 
density same as the blue cloud, 4 x 10 21 cm -2 , it can be estimated as ~1000M o with the 
radius of the 8pm ring of 1.7pc, the total mass of the accumulated layer is estimated to 
be ~2OOO-4OOOM 0 . This estimate is equivalent to ~ 50-95% of the mass of the ring 
component, 4.2 x 10 3 M o (see section 3.5), whereas the mass of the ionized gas is estimated 
as ~ 140 M 0 in section 3.5, which only accounts for less than 10% of the estimated mass 
of the compressed layer. These estimates suggest that ionization by the O star for the 
compressed layer created by the CCC is minor in mass, and that the remaining huge 
mass is now observable as the dense molecular ring surrounding the Hll region. This is 
consistent with our argument in the previous paragraph that the O star and the Hll region 
in RCW 120 are young, and suggests that the inner structures of the dense ring component 
were formed before the onset of the interaction with the Hll region. 
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By comparing the cavity depth 5 pc and the initial collision velocity of 30kms 1 , the 
lower limit of the timescale of the cavity creation and the subsequent O star formation in 


RCW120 can be estimated to be as short as ~0.2Myr, which is still consiste nt wit 
observational limit for the exciting O star in RCW 120, as discussed bv lMartins et al. 


i the 


(1201011 


that any age less than 5Myr is possible. This timescale may be larger like 0.4 Myr if we 
adopt 15kms _1 as the average relative velocity because the colliding velocity was probably 
decelerated during the collision as discussed above. We will investigate the relationship 


between the cavity depth and colliding velocity into more detail in the for 
(Torii et ah 2015b, in prep.). Considering the argument by Inoue & Fukui 


hcoming paper 


( 2013 1 that the 


effective Jeans mass and probably the mass of the forming stars are proportional to the 
third power of the colliding velocity, it is likely that the exciting O star in RCW 120 was 
formed in a short timescale at the beginning of the collision when the colliding velocity is 
still fast. The mass accretion rate to from the O star is estimated to be 0.5 x lO~ 4 M 0 yr _1 
for a stellar mass of ~2OM 0 by assuming that the accreting timescale is ~0.4Myr. Such a 
high mass accreti on rate is consistent with those revealed in the numerical calculations of 


CCC (llnoue & Fukui 


Tan fe McKee 


2003 


20131 ) and the theoretical studies on high-mass s tar formation (e.g., 


Krumholz et al. 


2009 


Hosokawa fe Omukai 


20091). A similar short 


timescale of ~ 0.1 Myr for the high-mass star formation is observationally derived in the 
super star clusters N GC 36 03 and West erlund 1 by careful optical studies of the young star 


by HST and VLT (Kudryavtseva, et al. 


studied with semi-analytical models by 


20121). T 


re short timescale in these clusters was 


Dib et. al. 


(120131 ). in which a strong increase in the 


star formation efficiency per unit time is required. 


We here summarize the CCC in RCW 120. About 0.4 Myr ago, the collision occurred 
between the large red cloud with size 12 pc and mass 5.1 x 1O 4 M 0 and the small blue 
cloud with size ~3.5pc and mass ~1OOOM 0 , which was accompanied by the V-shape cloud 
with 7.5 x 1O 3 M 0 . By assuming the decelerated colliding velocity of ~15kms -1 and the 
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cloud motion having an angle of 45° to the line-of-sight, the cavity was rapidly created in 
the red cloud in a short timescale of 0.4 Myr, forming the compressed layer with mass of 
2000-4000 Mq, whose number density increased to ~ 10 5 cm -3 from ~ 10 4 cm -3 through 
the collision. The cavity creation stopped inside the red cloud at the depth 5 pc, where the 
colliding small cloud was dissipated into the compressed layer. Massive clump formation 
and subsequent O star formation proceeded along with the cavity creation at a large mass 
accretion rate up to lO _4 M 0 yr _1 . After the birth of the O star, an Hii region with mass 
~200 M 0 was immediately created, and filled out the cavity and began interaction with 
the surrounding dense layer, which is now observed as the molecular ring. The dense 


comp ressed lay er was 


(IZavagno et al. 


20071 : 


porous, and so t 


Deharveng et al. 


re ionized gas partly ( ^25 %) streamed out though i t 


2009 


Anderson et al. 


2 014 


Walch et al. 


2011 


2012 ). 


The ionized gas just reached the opening of the cavity where the V-shape cloud is currently 
located, and now it is the very beginning of the erosion by the ionized stream which can be 
observed as the complementary distribution between 
assuming the velocity of the ionized gas of lOkrns" 1 (IZavagno et al. 


he blue cloud and t he Hii region. By 
20 o3 ), the timescale to 


fill the cavity with the ionized gas is calculated to be 0.5 Myr for the cavity’s depth 5 pc, 
which is consistent with our estimate of the timescale of the CCC in RCW120. 

The present CCC model in Figured^] predicts secondary CCC in RCW 120. As shown 
in Figure^c), the CO peak in the right side of the V-shape coincides with a local CO peak 
in the red cloud. In addition, broad velocity features which connect the red cloud and 
the blue cloud are found in the Hi emission (Figure[3]). The bridging feature in velocity 
m ay be u n ders t ood a s an intermediate velocity feature expected in an early stage of CCC. 


Haworth et al 


of CCC in 


(1201 51) d isc ussed that, by using the data of the numerical calculations 


Takahira et al. 


(2014j), the turbulent motion in the compressed layer can be 


observed in a position-velocity diagram as intermediate velocity features connecting two 
colliding clouds. On the other hand, one single broad velocity feature is identified with 
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the CO dataset at the rim of the V-shape cloud toward the direction of condensation 6 
with some YSOs (Figured^). These observational signatures may indicate the early stage 
of another collision between the V-shape cloud and the top of the red cloud as shown in 
Figurefl3l If it is correct, the coinciding CO peaks may correspond to the compressed layer 
by the collision. Unfortunately, our Mopra and ASTE observations did not cover these 
coinciding peaks and high velocity features, and in order to investigate the present idea and 
the possibility of the triggered star formation in future, further detailed observations of a 
large area are necessary. 


Formation of the second generation stars in the context of CCC is another issue 
to be addresse d. YSOs and their evolutional stages in RCW120 were investigated by 


Zavagno et al. 


(200 3) b ased on the Spitzer /GLIMPSE point source catalog. Subsequently, 


Deharveng et al 


(120091 ) made identification for the 138 sources whose positions and 24 pm 


brightnesses were iteratively measured from the Syitzer/M lPS 24 //in maps. The main 


difference between these two studies is that in 
distributed inside the 8 pm ring, while in 


Zavagno et al. 


Deharveng et al. 


( 20071 ) many YSOs are 


( 2009!) the re are almost no YSOs 


inside the ring. Using the carefully measured 24 pm stellar samples, 


D eharveng et al. 


(120091 ) 


applied various criteria for a wide range of wavelength from near-infrared to 24 pm, making 
their YSO identification firm and reliable, and we thus adopt the Deharveng’s catalog in 
the present study. It is, however, noted that a completeness issue possibly remains in the 
Deharveng’s catalog, particularly toward inside the 8 pm ring, because emission from warm 
dust grains heated by the O star is predominant there, making it difficult to make reliable 
identification of the YSOs (see Figured]). In fact, the statistical studies of the stellar 


distributions in t 

le Spitzer bubbles indicate s 

(ThomDson et al. 

2012 ; 

Kendrew et al. 

2012 ) 


20121). Analysing their numerical calculations of 


the triggered star formation by Hll region, 


Walch et al 


( 20131 ) pointed out that the YSO 


distributions are peaked around the ionization front, and discussed that the existence of 
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stars inside the bubbles seen in the observational studies can be explained as projection 
effects of three dimensional structures. Same as many other Spitzer bubbles, RCW120 
possibly harbors the YSOs inside the 8 /(in ring, although we have to await for future 
observations to better understand the YSO distribution in RCW 120. 


A most remarkable YSO feature identified in RCW 120 so far is seen in condensation 


1 , w. 


r ere a chain of YSOs ar e aligned parallel to the rim of the 8 /im ring (jZavagno et al 


200 ?]; 


Deharveng et al. 


20091 ). These authors discussed that the YSOs were formed via 


gravitational instabilities in the collected layer created by a pressure driven Hu region. 


IZavagno et al. 

2007; 

Deharveng et al. 

Oil 

o 

o 

CM 


2009|). On the other hand, another possibility is 


raised in the CCC scenario. As already discussed, the high-mass star was likely formed 
in a short timescale at the beginning of the collision when the colliding velocity was fast. 
Then, it is possible that the low-mass star formation may be triggered in a later stage when 
the colliding velocity becomes even smaller. It is interesting to note a speculation that 
the chain of YSOs located near the symmetry axis of the collision may have been formed 
under triggering in the final stage of the collision, where the cloud motion was significantly 
decelerated to a velocity much less than SOkrns" 1 . The low velocity will lead to smaller 
turbulent motion and smaller mass accretion rate in the compressed layer, providing a 
possible alternative to the conventional idea of triggering by the O star ionization. 


5. Summary 

The conclusions of the present work are summarized as follows; 

1. We performed CO ,7=1-0 and J= 3-2 observations toward the Galactic Hll region 
RCW 120, which has a beautiful infrared ring, with the NANTEN2, Mopra, and 
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ASTE telescopes, and have identified two molecular clouds associated with RCW120 
around radial velocities of —8kms _1 and —28kms _1 . The former (the red cloud) 
with total molecular mass of 5 x 10 4 M 0 has a ring-like component surrounding the 
Hu region in RCW 120. The latter (the blue cloud) with total molecular mass of 
8 x 10 3 M 0 has a V-shape at the north of the opening of the 8 //rn ring. Both of the 
red cloud and the blue cloud show high 12 CO <7=3-2/J=1-0 ratio of > 0.8 and high 
kinetic temperatures of > 30 K, suggesting that they are heated by the exciting O 
star and/or by the Hu region in RCW 120. These observational signatures suggest 
that, despite of a large velocity separation of 20kms _1 , the blue cloud as well as the 
red cloud is physically associated with the Hii region in RCW 120. 

2. The large velocity separation cannot be understood as a gravitationally bound system. 
An expanding Hu region is an idea to understand the origin of the molecular gas 
in RCW 120. In the assumption of uniform density distribution, Collect & Collapse 
(C&C) is expected to form a dense shell layer, which may explain the observed 
signature of RCW 120. Our observations, however, indicate no expanding motion 
in the velocity distribution of the molecular gas, which is inconsistent with what is 
expected in C&C. The theoretical studies on the non-uniform density distributions 
also do not provide a plausible explanation for the main observational signatures, i.e., 
the intensity distributions of the molecular ring in the red cloud, coexistence of the 
two clouds with a large velocity separation of 20kms _1 , and absence of the expanding 
motion. 


3. We postulate the scenario that a collision between the red cloud and the blue cloud 
triggered the formation of the exciting O star in RCW 120 in a short timescale of 


0.2-0.4Myr. The basic i dea of t 
was first given by 


ie_present cloud-cloud collision (CCC) scenario 


Habe fc Ohta (119921) . In this scenario, the velocity separation 
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of 20kms~ 1 can be interpreted as a projection of the colliding velocity. Once the 
collision occurred, a cavity was created in the red cloud, followed by the O star 
formation in the compressed layer at the interface of the two clouds. The strong UV 
radiation from the O star then illuminates the inner surface of the cavity, making it 
observable at 8 /irn as a ring structure. The present CCC scenario provides a possible 
explanation for the present observations, and if this is correct, the ROW 120 provides 
the second case of the CCC triggering the formation of one single O star next to M20. 
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A. CO velocity channel maps of the NANTEN2 and Mopra datasets 

CO <7=1-0 velocity channel maps of the NANTEN2 and Mopra datasets are presented 
in Figuresfl5l-fl9l at a velocity step of ~2.5kms~ 1 . In each figure, the 12 CO <7=1-0 emission 
obtained with NANTEN2 and Mopra are shown at the left panels and the middle panels, 
respectively, while 13 CO and C 18 0 J=l-0 obtained with Mopra are presented in the right 
panels in color and contours, respectively. 
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Fig. 1 .— A color composite image of RCW 120. Green, 


8 //in ( Beniamin et al. 


due, and rec 


2003). the Spitzer /MIPS 24 pm ( Carev et al. 


250 /mi ( Zavagno et al 


show the Spitzer/TRAC 


2009 ). and the Herschel 


2010) ■ Large cross indicates the exciting star, and filled red circles, 


filled white circles, and filled white triangles indicate the Class I. intermec 


Class II or flat spectrum, and Class II YSOs identified by 


Deharveng et al. 


iate Class I- 


(120091 ). Small 


crosses and labe 


Deharveng et al 


s indicate the cold dust condensations identified at 870 /mi observations by 


(2009!). Dashed lines show the observed region with Mopra. White contours 


show the outline of the 8 /mi ring which is used in the following figures to shows the outline 
of the 8 /mi ring, where the Spitzer 8 /mi image is median-filtered with a 9" x 9" window. 
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Fig. 2.— Large-scale 12 CO <7=1-0 distributions towards RCW120 are presented with the 
NANTEN2 dataset. The red cloud and The blue cloud are shown in (a) and (b), respectively. 
The observed area in the Mopra observations is here shown by white lines, (c) A comparison 
of the red cloud (image) and the blue cloud (white contours) is shown in a close-up view, 

in which the bridging feature at —23- 20kms _1 shown in Figure[3]is plotted in the thick 

black contours. The YSOs and dust condensations are plotted with the same manner as 
Figured] where the region which was used for the YSO identifications is shown bv black 
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Fig. 3.— Position-velocity maps for a very large area of RCW120. Contours show the 
averaged intensity of the NAN TEN2 12 CO data, and color image shows that of the SGPS 


Hi 21-cm data f McClure-Gri ffiths et ah 


200511 . Black solid line and dashed lines indicate 


position of the O star and the approximate extent of the 8 pm ring, respectively. The red 
arrow and blue arrows indicate the directions of the bridging features seen in CO and Hi, 
respectively. 
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Fig. 4.— 12 CO, 13 CO and C 18 0 ,7=1-0 distributions towards RCW120 with the Mopra 
dataset. The red cloud are shown in (a)-(c), while the blue cloud in (d)-(f). For the both 
clouds, 12 CO, 13 CO, C 18 0 J=l-0 emission are presented in the left, center, right panels, 
respectively. The unit of color bars is K km s -1 . Large cross indicates the O star, and 
the YSOs and dust condensations are plotted with the same manner as FigureHJ The gray 
contours in (a) are plotted at 90Kkms _1 and indicates the region which is used to estimate 
the molecular mass of the ring structure (see section 3.5). 
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Fig. 5.— (a-c) Comparisons of 13 CO distributions (color and contours) of the red cloud 
with (a) the Spitzer 8 fim emission (grayscale), (b) the SuperCOSMOS Her (grayscale), and 
(c) the 843 MHz SUMSS radio continuum emission (grayscale and thick blue contours). The 
integration range in velocity is from —12kms _1 to Okms -1 . (d-f) Comparisons of 12 CO 
distributions (color and contours) of the north component of the blue cloud in the same 
manner as (a)-(c). The integration range in velocity is from -^39kms -1 to —22kms _1 . 
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(a) Red cloud (-10 - +2 km/s) 
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Fig. 6.— i? 3 _ 2 /i—o maps of the red cloud (a) and the blue cloud (b) are shown in the color 
image with the 12 CO J=l-0 contours. The exciting stars, YSOs, and dust condensations 
are pl otted w ith the same manner as in Figure[fJ The names of the condensations given 


by 


Deharveng et ah 


(120091) are shown in (b). Regions A-F depicted with dashed circles 


are the regions that are used for the LVG analysis in Figure[7J (c, d) The region with high 
-R3-2/1-0 > 0.8 are shown in contours. The background is a composite image of the Herschel 
250 pm (red) and the Spitzer/IRAC 8pm (green). In (d), the outer boundary of the blue 
cloud is shown in yellow dashed contours. 
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Fig. 7.— LVG results for the regions A-F in Figure|6]are shown. i? 3 _ 2 /i_o and its error are 
shown with blue solid lines and filled blue area, while those of -R13/12 are shown in red. 
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Fig. 8.— Schematic images of the expanding shell model. Spatial distribution is shown in 
(a), while position-velocity maps expected in the model with different ratios of n ex p to Av 
are shown in (b) and (c). 
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Fig. 9.— (a) Longitude-velocity map of the Mopra 12 CO dataset. Vertical solid line and 
dotted lines indicate the position of the O star and approximate extent of the 8 /mi ring, (b) 
A close up of the red cloud in 13 CO (color) and C 18 0 (contours). 
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(a) first moment map _ 9 , 8 ™7 km/s (b) second moment map km/s 
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Fig. 10.— (a) The first moment map and (b) the second moment map of 13 CO ,7=1-0 in 
the red cloud. The contours shows the dense molecular ring in 13 CO J=l-0 integrated over 
—15-+2 kms -1 . 
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Fig. 11.— Azimuthal distributions of the red cloud with the NANTEN2 12 CO and 13 CO 
datasets are shown in (c), where the target regions of the plots are shown in (a) and (b) as 
a ring superimposed on the 12 CO and 13 CO integrated intensity maps. In (a) and (b), the 
inner boundary of the target ring is determined by fit to the 8 /ini image with a elliptical 
function, and the ring with 1.5 times larger size is then plotted as the outer boundary. The 
resulting major axes of the inner and outer rings are 1.6 pc and 2.4 pc, respectively, and the 
eccentricity is 0.24. The center of the ellipse is estimated to be (/, b ) = (348?257, 0?483), while 
the center of the azimuthal plots is taken at the position of the O star. In the azimuthal 
plots in (c), Az = 0° is defined at the point where the distance from the O star to the ring 
is smallest. 
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Fig. 12.— Schematics of the basic scenario of CCC introduced by 
section 4.2.2 for details). 
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Fig. 13.— Schematics of the evolution of the CCC in RCW120 are presented on the Z-Y 
plane in (a) and the X-Y plane in (b), where the X- and Y-axes are defined as the —90° -90° 
line and the 0°-180° line in Figures ITTT al and (b), and the Z-axis is along the line-of-sight. 
The origin of the coordinate system is taken at the exciting O star. 
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(a) Red cloud V L8R : -15 - +2 km s _1 


Black: 12 CO J=1-0, Red: ,3 CO J=1-0 
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Fig. 14.— CO integrated intensity distributions of the red cloud along the X-axis and the 
Y-axis are shown in (b) and (c), respectively. The black lines and the red lines respectively 
indicate the 12 CO and 13 CO. The NANTEN2 12 CO integrated intensity map is shown in (a), 
in which the X- and Y- axes are plotted with ticks. 
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Fig. 15.— Velocity channel maps of the CO emission at a velocity step of ~2.5kms~ 1 . The 
central O star is depicted by cross, and approximate boundary of the 8 pm ring is shown 
by white contours. The Mopra observing area is shown by white lines. (left column ) 12 CO 
,7=1-0 (color and black contours) obtained with NANTEN2. (center column ) 12 CO <7=1-0 
(color and black contours) obtained with Mopra. (right column) 13 CO <7=1-0 (color) and 
C 18 0 <7=1—0 (black contours) obtained with Mopra. The velocity ranges are shown at the 
left-top of the left panels in a unit of kms -1 . 
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Fig. 16.— Continued. 
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Fig. 17.— Continued. 
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Fig. 18.— Continued. 
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Fig. 19.— Continued. 











